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Abstract 

The pion scalar radius is given by (r|) = (6/n) J^ 2 ds 5s(s)/s 2 , 

with 5s the phase of the scalar form factor. Below KK threshold, 5s = 5 W , 
5-x being the isoscalar, S-wave nir phase shift. At high energy, s > 2 GeV , 5s 
is given by perturbative QCD. In between I argued, in a previous letter, that 
one can interpolate 5s ~ 5 n , because inelasticity is small, compared with the 
errors. This gives (r s ) = 0.75 ± 0.07 fm 2 . Recently, Ananthanarayan, Caprini, 
Colangelo, Gasser and Leutwyler (ACCGL) have claimed that this is incor- 
rect and one should have instead 5s — 5^ — it; then (r s ) = 0.61 ± 0.04 fm . 
Here I show that the ACCGL phase 5s is pathological in that it is discon- 
tinuous for small inelasticity, does not coincide with what perturbative QCD 
suggests at high energy, and only occurs because these authors take a value for 
5Ty(4m K ) different from what experiment indicates. If one uses the value for 
5 w (4m|f) favoured by experiment, the ensuing phase 5s is continuous, agrees 
with perturbative QCD expectations, and satisfies 5s — 5 n , thus confirming 
the correctness of my previous estimate, (r|) = 0.75 ± 0.07 fm 2 . 



Typeset with P H ys Ma TfrjX 



-THE SCALAR RADIUS OF THE PION- 



1. Introduction 



The quadratic scalar radius of the pion, (r|), is defined via the scalar form factor, Fs ;7r : 

F s ^(t) t ^ o FsAO){l + k(rs)t}, (1-1) 

where 

(tt( P )\ [m u uu(Q) + m d dd(0)] Hp')) = (27r)- 3 i^(t); (1.2) 

the (charged) pion states are normalized to (n(p)\iv(j/)) = 2p 5(p — p'), and t = (p — p') 2 . To one loop 
in chiral perturbation theory (ch.p.t.), (r 2 .) is related to the important coupling constant U by 

(rl) = J^{h-%}- (1-3) 

/tt ~ 93 MeV is the decay constant of the pion. 

An evaluation of (r§) was given some time ago by Donoghue, Gasser and Leutwyler;W we will 
refer to this paper as DGL. These authors found (we quote the improved result from the second paper 
in ref. 1) 

( r 2) DGL =0.61 ±0.04 fm 2 , r 4 = 4.4±0.2. (1.4) 

The error comes from experimental errors and the estimated higher order corrections. 
As noted in ref. 2, one can obtain the scalar radius from the sum rule 

d^, (1.5) 

AMI s 



= s / 

J 4 



where 8s(s) is the phase of Fs i7T (s), and is the charged pion mass. At low energy, 6s(s) = S n (s), 
where 6 n (s) is the phase shift for rnr scattering with isospin zero in the S wave. This equality holds 
with good accuracy up to the opening of the KK threshold, at s = &m 2 K ; for m K we take the average 
kaon mass, mx = 496 MeV. At high energy, s > 2 GeV 2 , one can use the asymptotic estimate that 
perturbative QCD indicates for 5s{s) (see below) and, between these two regions, what was considered 
in ref. 2 a reasonable interpolation, viz., Ss(s) ~ S n (s). One then finds, 

(r 2 ) =0.75 ±0.07 fm 2 , ^ = 5.4 ±0.5. (1.6) 

This is about 2 a above the DGL value, Eq. (1.4). 

The integral in (1.5) up to s = Am 2 K can be evaluated in a fairly unambiguous manner, and 
the contribution of the high energy region, s > 2 GeV 2 , although evaluated with different methods, is 
found similar in refs. 1, 2, 3. The conflictive contribution is that of the intermediate region, 



2GeV^ 



/ 

J4m 2 K 



t s 5 J^l. (1.7) 



In fact, very recently Ananthanarayan, Caprini, Colangelo, Gasser and Leutwyler,! 3 ! that we will denote 
by ACCGL, have challenged the result of ref. 2. Their main objection is that the Fermi-Watson final 
state interaction theorem does not guarantee that 5 n (s) and 5s(s) are equal, even if inelasticity is 
negligible; it only requires that they differ in an integral multiple of tt: 

6s(8)~6*(8) + Nir. (1.8) 

At 7T7T threshold, both 5 S and S n vanish, hence N = here. Below s = Am 2 K , continuity guarantees 
that the N in (1.8) still vanishes, as assumed in ref. 2. For 1.7 GeV 2 < s < 2 GeV 2 inelasticity is 
also compatible with zero. However, since this is separated from the low energy region by the region 
2itik < s 1 / 2 < 1.2 GeV, where inelasticity is not negligible, one can have 7^ 0. Actually, ACCGL 
conclude that 

6 s (s)~6 n (s)-7r, 1.1 GeV ~ s 1 ^ 2 ~ 1.42 GeV . (1.9) 
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According to ACCGL, this brings the value of (r§) back to the DGL number in (1.4). 

The remark of ACCGL leading to (1.8) is correct. Nevertheless, we will here show that their 
conclusion (1.9) is wrong. In fact, arguments of (A) Continuity of 5s(s) when the inelasticity goes 
to zero; (B) The value experiment indicates for the quantity 5- K (Am 2 K ); (C) The value SU(3) ch.p.t. 
implies for the real part of the KK scattering length and 8 n (Am 2 K ); and, finally, (D) The matching 
with the phase expected from perturbative QCD at high virtuality, all imply that the number N in 
(1.8) vanishes, therefore substantiating the claims of ref. 2. 

It should also be noted that the error analysis of DGL and ACCGL must be incomplete. With 
a correct error analysis, and even starting from their assumptions, DGL and ACCGL should have 
obtained a value for (r§) compatible with that in ref. 2, within errors. This is also discussed below. 

2. Some definitions 



Since we will only consider the S wave for isospin zero, we will omit isospin and angular momentum 
indices. We define a matrix for the partial wave amplitudes for the processes irn — > tttt, tttt — > KK(— 
KK -> tttt), and KK -► KK: 



f tttt — >K K fRK^KK 



( 



2i<5„ 



2i 



ve 2iS K _ l 

2i 



(2.1) 



Below KK threshold, the elasticity parameter is rj(s) = 1; above KK threshold 1 one has the bounds 
< r\ < 1 . We will also use a K-matrix representation of f : 



f={Q-V*K-Q-"S-i}-\ <,_(* °). 



(2.2) 



q a are the momenta, q± = \J s/A — M 2 , q 2 = \/s/4 — m 2 K . 
We may diagonalize f and find the eigenphases, 



gD = 



f = C{g c -i}- 1 C T , 
C = 



cot <S(+) 
cot<S(-> 



cos t) sin f 
- sin 9 cos ( 



(2.3a) 



We will define 6^ to be the eigenphase that matches 5^: S^^Am 2 ^) = 5 7r (4m| r ). Then, 



tan 8^ 



T ± VT 2 - 4A 



1 T + VT 2 - AA - 2 qi K n 



1/2 



(2.3b) 



L 2 +VT 2 - A A J 

T =q 1 K n + q 2 K 2 2, A = qi q 2 det K. 



This holds (near KK threshold) when K\\ > 0. For K\\ < 0, the (±) signs should be exchanged in 
the right hand side of the expression for tan 8^\ and the square roots in the expression for sin 6* get a 
minus sign. Near KK threshold, sin# and 8^(s) vanish with q 2 . If inelasticity were zero (77 = 1) the 
channels would decouple and one would have C = 1 and d (+ ^ = 8^, 5^ = 5k- 



In the present paper we will neglect coupling of tttt to states other than KK, for energies below 1.42 GeV. 



-2- 



-THE SCALAR RADIUS OF THE PION- 



The phase of the tttt — > rnr amplitude will play an important role in the subsequent discussion. 
We will actually use the phase 4> n defined by 



fn 



+ |/7r7r-»7T7r| ^ , < (j> n < 7T; 

I f 7T7T — >7T7T | , 7T < < 2VT. 



This definition has to be adopted to agree with the standard definition of the phase (shift) 8 for a 
purely elastic amplitude, given by / = sinJe" 5 , so that / = ±|/|e 1<5 with the (±) signs as for 4> n above. 
Using (2.1) one gets a simple expression for the tangent of (f) n : 

tan (f) w = 1 1 + i^-p (l + cot 2 8„) | tan 8 n . (2.4) 

For ease of reference, we also give here the expressions of phase shift and inelasticity in terms 
of the K-matrix: 



tan^ 



( gilggj detK + giKu 2 
- — — , s s 4%, 

1 + \q 2 \K 22 

2 9l [X 11+9 U 2 detK] {^ " ^ + ^^( detK ) 2 " 1 < 2 - fa > 
+ + q\Kl 2 + g ?g2 (detK) 2 + 1)2 _ i^^T | 5 g > 4m 2 K . 



(l + qiq2 detKr + ( qi K n - q2 K 22 r 2 

7? "V(l-^ 2 detK) 2 + ( (Zl K 11+ , 2 il' 22 ) 2 ' S ~ 4mK - (2 ' 5b) 

The connection with the scalar form factor of the pion comes about as follows. We form a 
vector F with Fs, n and the form factor of the kaon, Fs,k, and define the vector F' by 



F' = C T Q 1 / 2 F, F = (I s * V 

\*S,K J 

Then two-channel unitarity implies that 

F S ,n = q ^ /2 {(cose)\F{\e is(+> +(sin^)|^|e i5M }, 
Fs,k = q2 l/2 {(cos6)\F! 2 \e^ ) - {sm6)\F[\ e i5<+ > } . 

Near KK threshold or for small inelasticity, 8 s — 8^ ~ <f>^. 

3. The partial wave amplitudes from the experiment of Hyams et al. 



(2.6a) 



(2.6b) 



We will here follow DGL and ACCGL and take the partial wave amplitudes as measured by Hyams 
et al.,[ 4 l although later we will also discuss other sets of mr scattering data, as well as data' 5 ! on 
7T7T — ► KK. Hyams et alii give three representations for their data: an energy-independent phase shift 
analysis that yields the values of the phase shift 8 n (s), and of the elasticity parameter rj(s), from tttt 
threshold to s 1 ' 2 ~ 1.9 GeV; an energy-dependent parametrization of the K-matrix that interpolates 
these data in the whole range; and a second parametrization with a constant K-matrix that represents 
the data around KK threshold. 

For the second, Hyams et alii write [Eq. (12a) and Table 1 in ref. 4] 

K ab (s) = a a a fe /(s! - s) + (3 a (3 h j (s 2 - s) + j ab , (3.1) 
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where 

s }/ 2 =0.11 ± 0.15, s\ /2 = 1.19 ±0.01; 
ai =2.28 ±0.08, a 2 = 2.02 ±0.11; A = -1.00 ± 0.03, (3 2 = 0.47 ± 0.05; (3.2) 
7ii =2.86 ±0.15, 712 = 1.85 ±0.18, 722 = 1.00 ± 0.53. 

The numbers here are in the appropriate powers of GeV. 

In the energy range around KK threshold, 0.9 GeV < s 1/2 < 1.1 GeV, Hyams et alii (p. 148 
of ref. 4) find that their data may be represented by a constant K-matrix with 

if 11 = 1.0 ± 0.4 GeV -1 , K12 = -4.4 ± 0.3 GeV -1 , K 22 = -3.7 ± 0.4 GeV -1 . (3.3) 

The sign of K± 2 is undefined. We have chosen in (3.3) a sign opposite to that of Hyams et al.J 4 ! to 
agree with what the same authors get from the energy-dependent K-matrix; see below, Eq. (3.4). This 
is somewhat different from what (3.2) gives at KK threshold: evaluating K(s) with the central values 
in (3.2) one finds 

K xx (4m 2 K ) = -0.17 GeV -1 , K l2 (Am 2 K ) = -4.0 GeV -1 , K 22 (Am 2 K ) = -2.7 GeV -1 . (3.4) 

Before starting with the actual analyses it is perhaps convenient to remark that what follows 
from experiment is the energy-independent set of phase shifts and elasticity parameters. The energy- 
dependent representations are model dependent. This is particularly true of (3.1), where one makes 
the choice of a specific functional form; the results vary somewhat if using other parametrizations. 

4. The phase 



We will here consider the value of the phase 4> n (s) that follows from the experimental analysis given 
above. Although (j) n (s) is different from the quantities Ss(s) and <5( + )(s), which are the ones that 
intervene in the evaluation of the scalar form factor, they follow the same pattern. This was noted by 
ACCGL, who discuss 4> n in detail to illustrate their conclusions on 5s, and, indeed, it can be verified 
without too much trouble with the formulas of Sect. 2: explicitely for 5^ and to first order in q 2 or 
in e for 5 s (the exact result for the last requires solving two coupled integral equations). 

The advantage of <f> n is that it is given by the simple equation (2.4) in terms of the observable 
quantities 5 n , r\. This will allow us to simplify the discussion enormously; in particular, it will let us use 
simple parametrizations of 5 n , 77 above KK threshold, which is the region where there is disagreement 
between the evaluation of ref. 2 and DGL, ACCGL. This simplification is unnecessary in the sense that 
the results are almost identical to what one finds with the full K-matrix (that we will present later); 
but it permits us to exhibit, with great clarity, both the mechanisms at work and the issues involved. 

To calculate cb n around and above KK threshold we take 



5 w (s) = 7T + d(s) + -jj^c(s), 



V(s) =1 - e(s) 

and approximate, for 0.95 GeV < s 1 / 2 < 1.35 GeV, 



(4.1) 



d(s) — do — constant, c(s) = c = constant, 

2X M 2 -s „, , , ^ „ ( 4 - 2 ) 



*> = ( £l ^ + £ ^)^' M = 1-5 GeV . 
In the region immediately below KK threshold we replace q 2 by l^l in (4.1 
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FIGURE 1. The phases S w (s) (continuous line) and 4> w (s) (dashed line) cor- 
responding to a Type I solution. Note that, below KK threshold 5 w (s) = 
4>tt(s), hence fa(s) shows a very pronounced spike at s = 4m|. The asymp- 
totic phase (to be defined below) 5 as . is represented by the thick gray line. 



The energy-independent set of data in ref. 4 are well fitted with the numbers 2 

c = 5±l, ei = 6.4±0.5, e 2 = -16.8 ±1.6 (4.3) 

and we will leave the value of do (which is small) free for the moment. It will turn out that a key 
quantity in the analysis is the phase shift at KK threshold, J^^mf^ ), and we want to be able to vary 
this. 

4.1. The phase fa of DGL, ACCGL 

The authors of refs. 1, 3 take the K-matrix of Hyams et al.,1 4 ! with the central values as given in (3.2). 
What is important for us here is that this implies that the central value of 8^{Am 2 K ) is less than 180°: 

5^{Am 2 K ) = 175°. (4.4a) 

To reproduce this, we have to take do in (4.2) negative and equal to 

do = -0.087. (4.4b) 

Care has to be exercised when crossing the energy sq at which ^(s) equals n, which, with (4.3) and 
(4.4b), occurs at s^ 2 = 992.6 MeV, 

5 n (s = (992.6 MeV) 2 ) = tt, 

and where (2.4) is singular. For the moment, we will tackle this by starting below s and requiring 
continuity of the phase fa(s) across s . This we will call a solution of Type I, and is like what ACCGL 
find; indeed, the corresponding values of S n (s), fa(s), shown in Fig. 1, are practically identical to those 
in the Fig. 1 in ACCGL in the relevant region, around and above KK threshold. As can be seen in 
both figures, in the region s 1 / 2 ~ 1.35 GeV, where inelasticity is negligible, ^(s) and fa(s) differ by 
tt. 5s(s) and <5^ + ^(s) are very similar to fa(s) and thus also differ by tt from ^(s). 

2 We have actually followed the fit of ref. 6, which takes into account other data sets and is slightly below, both 
for 5^ and e, from what Hyams et al. give, at the upper energy range. 
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Figure 2. Fit to the / = 

and 1400 MeV, from ref. 6 



0, S'-wave inelasticity and phase shift between 950 



so that the formula used for <5. 



(0) 



5^ is slightly 



different from (4.1)], and data from refs. 4, 5, 8 and 11. The shaded bands cor- 
respond to 1 a variation in the parameters of the fits. The fit to the phase shift 
corresponds to d = 0. The difference between the determinations of n from 
nn — > tvtt (PY from data in the figure) and from tttx — > KK (PY alternative) 
is apparent here. 



This is the key remark of ACCGL: the phases 5- n {s) and ^(s), Ss(s) are not equal above 
s 1 / 2 ~ 1.1 GeV, but rather one has 

6s(s) - (s) ~ M») - <M S ) " tt ; s 1/2 > 1.1 GeV . 

This accounts for the difference between the results of refs. 1, 3 (DGL, ACCGL) and my previous 
results I 2 1 for the integral (1.7), hence for the different values of the scalar radius. 

The situation, however, is not as simple as ACCGL seem to believe. First of all, the inelasticity 
given in ref. 4 is much overestimated. After that paper was written, a number of experiments have 
appeared^ 5 ! in which the cross section 7T7t — ► KK was measured. Since there are no isospin-2 waves 
in 7T7T — > KK scattering, and the 7T7t — KK coupling is very weak for P, D waves, it follows that 
measurements of the differential cross section for 7T7t — > KK give directly 1 — n 2 with good accuracy. 
On the other hand, tttt scattering experiments like those of Hyams et alJ 4 l only measure the hit — > tttt 
cross section, so that rj is obtained less precisely here: not only the 7T7t cross section depends on both 
Jjr, 77, but other waves (notably S2, P and DO) interfere. Thus, these more recent, 7T7t — > KK based, 
experimental values^ 5 ! for rj are much more reliable than the older ones, in particular than those of 
ref. 4. 

The value of the inelasticity the experiments in ref. 5 give is about a third of what (4.3) 
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1.0 1.1 1.2 1.3 1.4 GeV 



FIGURE 3. The phases S w (s) (continuous line) and </>tt(s) (dashed line) cor- 
responding to a solution of Type Id. The spike that appeared in Fig. 1 is now 
accompanied by a jump of <^,r(s), at s = so- 



indicates: r\ can be fitted withJ 6 ! 

d = 2.4 ±0.2, e 2 = -5.5±0.8. (4.5) 

The difference is shown graphically in Fig. 2. 

If we now use (4.5) instead of (4.3) to calculate <f>w(s), keeping ^(s) fixed, a surprising result 
occurs: 4> n (s) does not become closer to S n (s) above the point so; on the contrary, it moves closer to 
6 W — 7r. In fact, one can decrease the inelasticity to zero, e(s) — > 0, keeping 5 n (s) fixed, and one finds 
that 

M*)-<M*), s 1 ' 2 < 4 /2 = 992.6 GeV; 
w (s)_,^( s )_ 7rj s 1 / 2 > s 1 q 2 = 992.6 GeV . 

That is to say: contrary to physical expectations, the limit of zero inelasticity does not coincide with 
inelasticity zero for, if we set e(s) = 0, then 5 7T {s) and ^{s) should be identical. This phenomenon 
was noticed by ACCGL who, however, failed to attach to it the due importance. As a matter of fact, 
the situation is even more complicated, as will be shown below: if we leave r\ fixed but vary do in (4.2) 
across zero to a positive number, however small, the resulting cf> n is not continuous when d crosses 
zero: it jumps by it. 

What is the reason for this peculiar behaviour of 4>tP- It is not difficult to identify: Eq. (2.4) 
does not determine 4> w , but only its tangent. Thus, <j>^ is only fixed up to a multiple Nir. N may be set 
to zero below the point sq where 5- K {s) crosses tt, by requiring that <?!) 7r (4m| ( -) = 5 v {Am K ) and continuity 
above this. However, Eq. (2.4) shows that tan (j) n (s) is discontinuous when s crosses so- Therefore, we 
may well add n to the (j>^{s) of the Type I solution found above, in the region s > s , since this does 
not change its tangent. We then find what we call a solution of Type Id (d for discontinuous), depicted 
in Fig. 3. In this case, ^(s) is not continuous across s , but does tend 3 to ^(s), for all values of s, 
when the inelasticity tends to zero. It is also continuous (in the mean) for do around zero. ACCGL 
appear to be unaware of the existence of solutions of Type Id. 

It is not clear which of the two solutions, Type I or Type Id, should be considered correct: 
both types look awry. In fact, we will show that both Type I and Type Id are, with all probability, 

3 To be precise, one should remark that this limit applies in the mean; the isolated point so remains singular. 
Convergence in the mean, however, is sufficient to ensure convergence of integrals involving <j> w . 
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FIGURE 4. The phases S w (s) (continuous line) and 4> w (s) (dashed line) cor- 
responding to the solution of Type II. As in Fig. 1, the thick gray line is the 
asymptotic phase (5 as .- 



spureous solutions, artifacts due to the use of the parametrization (3.1), (3.2) over too wide a range, 
and with too little experimental information. 

4.2. The correct <fi n 

We next repeat the calculations of the previous section, but we will now assume that 5 3r (4m| f ) is larger 
than 7T, so that do is positive. To get this it is sufficient to alter a little the parameters in (3.2). For 
example, if we move only one parameter by 1 c, just replacing in (3.2) 

ai -^ 2.20 = 2.28 - 0.08, (4.7) 

then 5 n (4m 2 K ) becomes 185°. Note that ^(s) is almost unchanged by this, as may be seen by comparing 
Figs. 1 and 4. The only important effect of the change in (4.7) is to push 8^{Am 2 K ) from a bit below to 
a bit above 180°; but then, this is a key point, as we will see. 

A value for 8 n (4m 2 K ) above 180° follows also for si = 0, 711 = 3.0 (as in ref. 7), which values 
are both less than 1 a off the central values in (3.2). In fact, a value S^Arn^) > 180° can already be 
obtained with only a \ a change, 

ai -> 2.24 = 2.28 - 0.04. 

Thus, a value S 7T (4:m'j <: ) > 180° is perfectly compatible with the energy-dependent parametrization of 
Hyams et alii, Eqs. (3.1), (3.2), when errors are taken into account. 

We will use (4.7) for simplicity in the discussion and will thus repeat the calculations with 

S w (4m 2 K ) = 185° , d = +0.087. (4.8) 

In the present case, and as is obvious from (2.4), </» 7r (s) is never singular and it stays above 5 n (s), up 
to the energy s 1 / 2 ~ 1.3 GeV where 5 n (s) crosses 37r/2, remaining close to it afterwards. 4 

This property is actually quite general, not tied to the specific approximations (4.2), (4.8), 
and depends only on the fact that S n (s) is an increasing function of s and that <5 7r (4m 2 ^) > it. This is 
all we need for 4> n . To get the analogous property for , ds we also require that detK(4m|-) < 0, 

4 In fact, over the whole range, the difference between (f)^ and 5^ is smaller than the experimental errors of the 
last: compare Figs. 2 and 4. 
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Figure 5. The phases 

S-k(s) (continuous lines) 
and <5( + )(s) (dashed lines) 
evaluated in the K-matrix 
formalism, Eq. (3.1), 
with the central values 
of the parameters given 
in (3.2) (Type I) or (4.7), 
Type II. [The two lines 
for 5„ correspond also to 
(3.2), (4.7)]. The asymp- 
totic phase 5as. (thick 
gray line) is also shown. 



something that is amply satisfied with the parameters of (3.2), (3.3) or, more generally, if, as implied 
by SU(3) ch.p.t., one has taxiSx < near KK threshold (see below). 

A set of phases with these properties we will call a solution of Type II. In the specific case 
(4.8) we find the S n (s), (p n (s) depicted in Fig 4. Note that 5 n (s) and (?V( S ) are near each other all the 
time, as one expects physically since the inelasticity is small; this is particularly important in view of 
the results of refs. 5. Unlike what happened in solutions of Type I, or Type Id, the phase 4> w (s) is now 
a smooth function both of e(s) and of s. 

These results are not new. They were amply discussed more than thirty years ago, in connection 
with the eigenphases <$W(s), by the present author in ref. 7. There it was noted that, by going from the 
values of the K-matrix parameters in (3.2) to values like those in (4.7), the eigenphase <$(+)(s) changes 
from a fast decrease above the KK threshold, diverging from <J T (s) by ~ tt (as does (f> n in a solution 
of Type I, see Fig. 1), to increasing above KK threshold with increasing s, staying close, but a bit 
above, <5„-(s) (again, as does t/v in a solution of Type II, Fig. 4). The reader may compare our Figs. 1, 4 
here with Fig. 2 in ref. 7. In ref. 7 the M-matrix (M = K _1 ) parametrization of experimental data of 
Protopopescu et alJ 8 l is also considered, and the same phenomenon is observed (Fig. 1 in ref. 7). 

We give in Fig. 5 the eigenphases 5^ corresponding to Type I and Type II solutions. Here 5 V , 
as well as the eigenphase d^ + \ are evaluated with the K-matrix formalism, Eq. (3.1). For Type I we 
took the parameters (3.2); for Type II, those in (4.7). Our Fig. 5 here agrees with the corresponding 
parts of Figs. 1, 2 in ref. 7. 

5. The value of S^Arn^) 

As is obvious from the previous discussion, a key quantity in this analysis is the tttt phase at KK 
threshold, 5 7r (4m 2 < ). If this is smaller than tt, we have a situation of Type I; if, on the contrary, 
8^{4rn K ) > tt, we have a solution of Type II and, in particular, we can approximate 



as was done in ref. 2. 

It should be clear that the parametrization (3.1), (3.2) is not a good guide to find the value of 
5 7T (Am 2 K ). Not only d^^m 2 ^) crosses 180° when varying the parameters in (3.2) within their errors (as 
we have shown before) but, more to the point, (3.1) was devised to furnish an approximate represen- 
tation of <5,r(s), rj(s) in the whole range AM 2 to 1.9 2 GeV 2 . This may easily create local distortions; 
and, in fact, such distortions are expected. The inelasticity of ref. 4 is overestimated, as proven by 



6 s (s) ~ 6( + \s) ~ ^(s)^6 n (s), 
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the more precise measurements of ref. 5: this will influence the phase 8^ above f GeV, hence, via the 
parametrization, around KK threshold. Such a distortion also occurs in the evaluation of Au et al.,' 9 ' 
who make a fit to rj and d n , based on data of ref. 4, over the whole energy range, which fit leads to a 
value of 5- K {Am 2 K ) smaller than 180°: see Fig. 4 in ref. 9. We certainly need something more precise in 
the vicinity of the KK threshold, since 5 vr (4m|-) is so near ir. 

For this we have several possibilities: the constant K-matrix fit around KK threshold of Hyams 
et al.;I 4 l the energy-independent analysis of this same reference; the results of other experiments; or 
certain theoretical arguments. As for the first, if we take the values K ab in (3.3), obtained from a fit 
to data from 0.9 GeV to 1.1 GeV, we find 

5n(Am 2 K ) = 205 ±8°, (5.1) 

3 a above 180°. A value above 180° is, of course, also found if interpolating the energy-independent 
analysis of ref. 4. The data of Protopopescu et alJ 8 l are not sufficiently precise to discriminate whether 
S^iim 2 ^) is below or above 180°: for some of the solutions in ref. 8, S^^m 2 ^) is below, and for others 
above 180°, but in all cases, the errors cover the value 180° (at any rate, the elasticity parameter of 
Protopopescu et al. is incompatible with nir — > KK results). However, a value clearly above 180° is 
found if extrapolating downward the experimental results of ref. 10 (the phase shift is only measured 
for s 1 / 2 > 1 GeV). This gives 5 S n (Am 2 K ) = 203 ± 7° , including estimated systematic errors. A value 
5 7r (4m^ r ) > 180° is also found in all five solutions of Grayer et al.:' 11 ! cf. Fig. 31 there. Finally, Kamihski 
et alJ 11 ! find 5 n (4:m 2 x ) = 190 ± 25° . The experimental information thus very strongly favours a value 
S 7T (4m 2 < ) > 180° , and hence a solution of Type II. 

There are two other independent, theoretical arguments in favour of a solution of Type II. The 
first is based on chiral SU(3) calculations: unitarized SU(3) ch.p.t. produces central values of 5 n (4m 2 K ) 
above 190° (with a value around 200° favoured; see for example ref. 12). Moreover, in Type II solutions, 
with the parameters in (4.7), one has a real part of the KK scattering length a r (KK) ~ — 0.46M" 1 , 
in agreement with the unitarized current algebra (ch.p.t.) result that gives a r (KK) ~ — 0.5M" 1 . 

The second, more serious indication, is that the phase Ss(s) — (finis) for Type II solutions joins 
smoothly the result furnished by the perturbative QCD evaluation of 5s(s), while a Type I solution 
5s(s) lies clearly below. We now turn to this. 

6. The phase 5s(t) at large t from QCD 



Using the evaluations in ref. 13 it is easy to get that, to leading order in the QCD coupling a s , one has 

47r[rr4(^)+m^ 2 )]C F a>\ (61) 

where, neglecting quark and pion masses, 

Here v 2 is the renormalization point and & is the partonic wave function of the pion, defined by 

(2tt) 3 / 2 (0|5 : J(0)7A75^ • ■ • D^O) : \tt( P )) = i n+ Ww • • -p„„ C (£, A 

Jo 



5 For the favoured solution in ref. 10 which, incidentally, is the one with values of r?(s) more compatible with 
measurements based on rnr —> KK. For other solutions S 7T (4m K ) is even larger, except for one that yields a 
value near 180° . This last one, hovever, has an elasticity parameter incompatible with mv — > KK results. 



FsAt) 
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The Dp are covariant derivatives, and S means symmetrization. The function & is the same that 
appears in the evaluation of the vector form factor, and thus^ 13 ! 

m» 2 ) - £(1-66^2/.- (6.3) 

If we input (6.3) into (6.2) we get a divergent result. This divergence may be traced to the fact that we 
have neglected quark and pion masses, and may be cured by defining the form factor not for external 
momenta p 2 = p' 2 = 0, but with p 2 = p' 2 = to, to being a fixed number; for example, we could take 
t = M 2 . Then we choose v 2 = — t (for spacelike t) and find the asymptotic behaviour 

FsAt) „ ^n[ml(-t)+m 2 (-t)]C F f 2 a s (-t)log(-t/t ) C[ mj (-t) + m 2 (-t)]f 2 

' t^oo — t —t 

with C = 5767r 2 CV/(33 — 2rij), and rif is the number of quark flavours, that we take equal to three. 

Unfortunately, the value of the constant C is changed when higher order corrections are in- 
cluded. These have the same structure as (6.4), with higher powers [a s (— t) log(— t/to)] n which are 
not suppressed at large t. Therefore, the constant C gets contributions from all orders of perturbation 
theory with the result that the final value is unknown. However, it is very likely that the structure 
[(Constant) x J2m 2 (— t)/t] remains. This is sufficient to get a prediction for the asymptotic phase: 

feW = ' j 1 + to^5)}. = 12/(33 -2n,). (6.5) 

Here A is the QCD parameter; in our calculations here we have allowed it to vary in the range 
0.1 GeV 2 < A 2 < 0.35 GeV 2 . 5 as .(s) is the phase plotted in Figs. 1, 4, 5 where it is seen very clearly 
that it is consistent with Type II solutions, but not with the Type I solution of ACCGL. 

7. Conclusions 

There are other methods for finding directly I4, of which we only mention two. One can evaluate on 
the lattice the dependence of the quark condensate on the quark masses; I 14 ! or one can fit I4 to the 
experimental tttt scattering lengths and effective range parameters obtained from experimental data,! 6 ! 
using ch.p.t. to one loopJ 15 ! The results are summarized below, where we also repeat the results of 
refs. 1, 2, 3: 

4.4 ±0.2 [refs. 1, 3] 
( 5.4 ±0.5 [ref. 2] 

' 4 ~^ 4.0 ±0.6 [lattice calculation, ref. 14] [ -'- L) 
7.2 ± 0.7 [fitting a\ I] , b { Q I] ', b ly ref. 15]. 

This is inconclusive; lattice calculations are known to suffer from large systematic errors, and the 
number following from the fit to experimental data is affected by higher order corrections, which the 
evaluation in ref. 15 does not take into account. We have to fall back on our previous discussion, 
involving the phase of the scalar form factor. 

In this case, and as we have shown, we have two types of solution: Type I, that occurs when 
5 7T (4 : m 2 < ) < 7r, and Type II, when S 7T (4:m 2 < ) > n. The correctness of a solution of Type I, which is 
the one used in the evaluations of DGL, ACCGL is very unlikely: the experimental indications^ 4 ' 10,11 ! 
favour values 5 v {Am 2 K ) > it. Moreover, in Type I solutions one has a discontinuous phase when 
the inelasticity tends to zero. Type I solutions also exhibit a phase 4> w which is not continuous when 
d^^m 2 ^) moves around n. Finally, Type I solutions give a phase Ss{s) rather different from what 
perturbative QCD suggests, Eq. (6.5), at large s. We think that Type I solutions are spureous, un- 
physical solutions, which appear only because one tries to fit, with too simple a formula, and without 
enough experimental information, the whole energy range from tttt threshold to 1.9 GeV, which dis- 
torts the results in the region of KK threshold. This last conjecture is confirmed by the evaluations of 
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MoussallamJ 16 ! This author uses, like DGL, ACCGL, fits that represent the quantities 8^ and 77 over the 
whole energy range; in particular, the fit of Au et alJ 9 l Such parametrization gives 8 1T {Am 2 K ) ~ 173° , 
hence a Type I solution and thus, not surprisingly, Moussallam finds a value for (r|) similar to that of 
DGL. 

Although this is not very important, because the very starting point of DGL, ACCGL (a Type I 
solution) is unlikely to be correct, one may question the methods of error analysis of these authors. As 
we discussed above, a value 5^{^m 2 K ) > 180° is obtained if replacing a\ — > 2.28 — 0.04, i.e., moving 
only I a off the central value in the fits of Hyams et al.,l 4 l Eq. (3.2) here. Variation within errors of 
their parameters should have taken DGL, ACCGL to a Type II solution and, therefore, their error for 
(r|) should have comprised the value found with a Type II solution. With a complete error analysis 
DGL, ACCGL should have got 6 (rf) = 0.61^;^ fm 2 . 

For a Type II solution, on the other hand, the value of 8 v (4m%-) > n agrees with what 
experiment indicates; the phases ^(s), 8^(s) and 8s (s) are continuous both in s and when the 
inelasticity goes to zero; and the phase 5s(s) agrees well with what perturbative QCD suggests at 
large s. We conclude that a situation of Type II is by far the more likely to be correct, thus confirming 
the validity of the approximations in ref. 2; in particular, the estimate 

(rf) = 0.75 ±0.07 fm 2 . (7.2) 

A last question is whether one can improve on the evaluation in ref. 2. This is very unlikely, 
for the contribution of the region Am 2 K < s < 2 GeV, Eq. (1.7). First of all, the incompatibility of the 
central values for r\ in analyses based on tt-k — > irir scattering^ 4 ' 10,11 ! with what one finds in tttt — > KK 
experiments, I 5 ' implies that the phase 8 n obtained from tttt — > tttt scattering must be biased. And, 
secondly, to find the eigenphases 8^ and mixing angle 9 which are necessary to disentangle the 
form factors Fs^, Fs,k [cf. Eq. (2.6)], one requires, as discussed in detail in ref. 7, experimental 
measurements of the three reactions rnr — > tttt, tttt — > KK, KK — > KK. Failing this, we are only left 
with approximate evaluations, like those in ref. 2. 



6 Note that the converse is not true, in the sense that we do not have to enlarge the errors to encompass the 
DGL number: while it is true that the parametrization (3.1), (3.2) is compatible with both a solution of Type I 
and one of Type II, we have shown in Sect. 5 that the experimental data point clearly to S 7 ,(4mj ( ) > 180° , 
hence a solution of Type II, that SU(3) ch.p.t. calculations also indicate a solution of Type II and, finally, in 
Sect. 6, we have argued that only a solution of Type II is compatible with the asymptotic behaviour indicated 
by perturbative QCD. 
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